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ABSTRACT 

Context. Earth-sized planets around nearby stars are being detected for the first time by ground-based radial velocity and space-based 
transit surveys. This milestone is opening the path toward the definition of missions able to directly detect the light from these plan- 
ets, with the identification of bio-signatures as one of the main objectives. In that respect, both the European Space Agency (ESA) 
and the National Aeronautics and Space Administration (NASA) have identified nulling interferometry as one of the most promising 
techniques. The ability to study distant planets will however depend on the amount of exozodiacal dust in the habitable zone of the 
target stars. 

Aims. We assess the impact of exozodiacal clouds on the performance of an infrared nulling interferometer in the Emma X-array 
configuration. The first part of the study is dedicated to the effect of the disc brightness on the number of targets that can be surveyed 
and studied by spectroscopy during the mission lifetime. In the second part, we address the impact of asymmetric structures in the 
discs such as clumps and offset which can potentially mimic the planetary signal. 

Methods. We use the DarwinSIM software which was designed and validated to study the performance of space-based nulling in- 
terferometers. The software has been adapted to handle images of exozodiacal discs and to compute the corresponding demodulated 
signal. 

Results. For the nominal mission architecture with 2-m aperture telescopes, centrally symmetric exozodiacal dust discs about 100 
times denser than the solar zodiacal cloud can be tolerated in order to survey at least 150 targets during the mission lifetime. 
Considering modeled resonant structures created by an Earth-like planet orbiting at 1 AU around a Sun-like star, we show that this 
tolerable dust density goes down to about 15 times the solar zodiacal density for face-on systems and decreases with the disc inclina- 
tion. 

Conclusions. Whereas the disc brightness only affects the integration time, the presence of clumps or offset are more problematic and 
can seriously hamper the planet detection. The upper limits on the tolerable exozodiacal dust density derived in this paper must be 
considered as rather pessimistic, but still give a realistic estimation of the typical sensitivity that we will need to reach on exozodiacal 
discs in order to prepare the scientific programme of future Earth-like planet characterisation missions. 

Key words. Instrumentation: high angular resolution - techniques: interferometric - circumstellar matter 



1. Introduction 

The possibility of identifying habitable worlds and even biosig- 
natures from extrasolar planets currently contributes to the 
growing interest about their nature and properties. Since the 
first planet discovered a round another solar- type star in 1995 
(May or & Ouelozl [l995). nearly 370 extrasolar planets have 
been detected and many more are expected to be unveiled by 
ongoing or future search programmes. Most extrasolar planets 
detected so far have been identified from the ground by indi- 
rect techniques, which rely on observable effects induced by 
the planet on its parent star. From the ground, radial velocity 
measurements are currently limited to the detection of planets 
about 2 times as massive as the Earth in o rbits around Sun- 
like and low-mass stars dMavor et al.l 120091) while the transi t 
method is limited to Neptune-sized planets dGillon et al.ll2.Q07b . 
Thanks to the very high precision photometry enabled by the 
stable space environment, the first space-based dedicated mis- 
sions (namely CoRoT and Kepler) are now expected to reveal 
Earth-sized extrasolar planets by transit measurements as well. 
Launched in 2006, CoRoT (Convection Rotation and planetary 



Trans i ts) has detected its fi r st extrasolar planets (e.g jBarge et al.l 
2008; lAlonso et ail 120081; iLeger et all IsubmJ) and is expected 
to unveil about 100 transiting planets down to a size of 2R ffi 
around GOV stars and 1 . 1 R Q ar ound MOV stars ov er its entire 
lifetime for short orbital periods (Mout ou et alj 2005). Launched 
in 2009, Kepler will extend the survey to Earth-sized planets lo- 
cated in the h abitable zone of about 10 5 main sequence stars 
dBorucki et al.l 12007b . After 4 years, Kepler should have dis- 
covered several hundred of terrestrial planets with periods be- 
tween one day and 400 days. After this initial reconnaissance 
by CoRot and Kepler, the Space Interferometry Mission (SIM 
PlanetQuest) might provide unambiguously the mass of Earth- 
sized extrasolar planets orbiting in the habitable zone of nearby 
stars by precise astrometric measurements. With CoRoT and 
Kepler, we will have a large census of Earth-sized extrasolar 
planets and their occurrence rate as a function of various stel- 
lar properties. However, even though the composition of the up- 
per atmosphere of trans iting extrasolar planets can be probed 
in favorable cases (e.g., Richards on et al.ll2007l) . none of these 
missions will directly detect the photons emitted by the planets 
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which are required to study the planet atmospheres and eventu- 
ally reveal the signature of biological activity. 

Detecting the light from an Earth-like extrasolar planet is 
very challenging due to the high contrast (~10 7 in the mid- 
IR, ~10 10 in the visible) and the small angular separation (~ 
0.5 /irad for an Earth-Sun system located at 10 pc) between the 
planet and its host star. A technique that has been p roposed to 
overc ome these difficulties is nulling interferometry (Bra cewelfl 
1 1978b . The basic principle is to combine the beams coming from 
two telescopes in phase opposition so that a dark fringe appears 
on the line of sight, which strongly reduces the stellar emission. 
Considering the two-telescope interferometer initially proposed 
by Bracewell, the response on the plane of the sky is a series 
of sinusoidal fringes, with angular spacing of A/b. By adjusting 
the baseline length (b) and orientation, the transmission of the 
off-axis planetary companion can then be maximised. However, 
even when the stellar emission is sufficiently reduced, it is gen- 
erally not possible to detect Earth-like planets with a static array 
configuration, because their emission is dominated by the ther- 
mal contribution of warm dust in our solar system as well as 
around the target stars (exozodiacal cloud). This is the reason 
why Bracewell proposed to rotate the interferometer so that the 
planetary signal is modulated by alternatively crossing high and 
low transmission regions, while the stellar signal and the back- 
ground emission remain constant. The planetary signal can then 
be retrieved by synchronous demodulation. However, a rapid ro- 
tation of the array would be difficult to implement and as a re- 
sult the detection is highly vulnerable to low frequency drifts in 
the stray light, thermal emission, and detector gain. A number 
of interferometer configurations with more than two collectors 
have then been proposed to perform faster modu lation and over- 
come this problem by using phase chopping dAngel & W oolf 
119971: iMennesson & Mariottilll997t lAbsill 120011) . The principle 
of phase chopping is to synthetize two different transmission 
maps with the same telescope array, by applying different phase 
shifts in the beam combination process. By differencing two dif- 
ferent transmission maps, it is possible to isolate the planetary 
signal from the contributions of the star, local zodiacal cloud, 
exozodiacal cloud, stray light, thermal, or detector gain. Phase 
chopping can be imple mented in var ious ways (e.g. inherent 
and internal modulation, lAbsill 120061) . and are now an essen- 
tial part of future space-based life-fi nding nulling interfer ome- 
try missions such as ESA's Darwin (Fridl und et al|l2006l) and 
NASA's Terrestrial Planet Finder (TPF. lLawson et alj20 08). The 
purpose of this paper is to assess the impact of exozodiacal dust 
discs on the performance of these missions. After describing the 
nominal performance of Darwin/TPF, the first part of the study 
is dedicated to centrally symmetric exozodiacal discs which are 
suppressed by phase chopping and only contribute through their 
shot noise. If they are too bright, exozodiacal discs can drive the 
integration time and we investigate the corresponding impact on 
the number of planets that can be surveyed during the mission 
lifetime. In the second part, we address the impact of asymmet- 
ric structures in the discs (such as clumps and offset) which are 
not canceled by phase chopping and can seriously hamper the 
planet detection process. 

2. Darwin/TPF overview 

Considerable effort have been expended in the past decade by 
both ESA and NASA to design a mission that provides the re- 
quired scientific performance while minimizing cost and techni- 
cal risks. After the investigation of several interferometer archi- 
tectures, these efforts culminated in 2005-2006 with two parallel 
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Fig. 1. Representation of the Darwin/TPF space interferometer 
in its baseline "Emma X-array" configuration (Lege r~& Herbstl 
120071) . It includes 4 telescopes and a beam combiner spacecraft, 
deployed and observing at the Sun-Earth Lagrange point L2. 
At any given time, it can observe an annular region on the sky 
between 46° and 83° from the anti-solar direction. During one 
Earth year, this annulus rotates and gives access to almost all 
regions of the celestial sphere. 



assessment studies of the Darwin mission, carried out by EADS 
Astrium and Alcatel-Alenia Space. Two array architectures have 
been thoroughly investigated during these industrial studies: the 
four-telescope X-arr ay and the Three-Telescope Nuller (TTN 
iKarlsson et al.ll2004b . These studies included the launch require- 
ments, payload spacecraft, and the ground segment during which 
the actual mission science would be executed. Almost simulta- 
neously, NASA/JPL initiated a simi lar study for the Te rrestrial 
Planet Finder Interferometer (TPF-I. ILawson et"aI1l2008l) . These 
efforts on both sides of the Atlantic have finally resulted in a con- 
vergence and consensus on mission architecture, the so-called 
non-coplanar or Emma-type X-array (represented in Fig.[TJ. 

2.1. Instrumental concept 

The baseline design consists of four 2-m aperture collector 
spacecraft, flying in rectangular formation and feeding light to 
the beam combiner spacecraft located approximately 1200 m 
above the array. This arrangement makes available baselines up 
to 70 m for nulling measurements and up to 400 m for the gen- 
eral astrophysics programme (constructive imaging). Note that 
the size of the collecting apertures has not yet been fixed and 
will influence the final cost of the mission. The optical layout 
separates the nulling and imaging functions, the shortest base- 
lines being used for nulling and the longest ones for imaging. 
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Table 1. Instrumental parameters considered in this study for 
Darwin/TPF. 



Telescope 1 



Instrumental parameters Emma X- Array Design 



Max. baselines [m] 
Telescope diameter [m] 
Field of regard 
Optics temperature [K] 
Detector temperature [K] 
Quantum efficiency 
Instrument throughput" 
Science waveband [yum] 
Modal filtering [pim] h 
Spectral resolution 
Instrumental stability^ 

- rms OPD error [nm] 

- rms amplitude error 



400 x 67 
2.0 

46° to 83° 
40 
8 

70% 
10% 
6.0-20.0 
6.0-11.5/11.5-20.0 
60 

1.5 
0.05% 



8 excluding ideal beam combiner losses and coupling efficiency. 
b see section [3721 
c see section |3~4"1 



This configuration has the advantage of allowing optimal tuning 
of the shorter dimension of the array for starlight suppression 
while keeping a significantly longer dimension to provide a rapid 
modulation of the planet signal as the array rotates. The X-array 
design is also appropriate to implement various techniques for 
removing instability noise, which is one of the dominant noise 
contributor (see appendix©. The assessment studies settled on 
an imaging to nulling baseline ratio of 3:1, based on scientific 
and instrument design constraints. A larger ratio of 6:1 could 
nonetheless improve performance by simplifyi ng noise r educ- 
tion in the post-processing of science images dLavl 12006). The 
optical system architecture is represented by the block diagram 
in Fig. [2] with the following elements in the optical path: 

- four spherical primary mirrors located on a virtual 
paraboloid and focusing the beam at the paraboloid's fo- 
cal point. The virtual parabola focal length has been set to 
1200 m, resulting from a trade-off between differential po- 
larisation effects and inter spacecraft metrology capability, 
as well as to enable the implementation of longer baselines 
for an imaging mode; 

- transfer optics consist of all the equipment needed to col- 
lect and redirect the incoming beams towards fixed direc- 
tions whatever the interferometer configuration. This in- 
cludes mainly tip-tilt mirrors to handle array reconfiguration, 
a derotator to handle the array rotation and a common two or 
three mirror telescope to achieve beam collimation; 

- optical delay lines (ODL) to adjust the optical path differ- 
ences (OPD); 

- fast steering mirrors to correct for tip-tilt errors; 

- deformable mirrors to compensate for quasi-static errors 
such as defocus and astigmatism which are due to the use 
of spherical mirrors; 

- achromatic n phase shifters to produce the destructive inter- 
ference of on-axis stellar li ght (using Fresnel Rhombs for 
instance, Mawet et al. 20071); 

- dichroic beam splitters to separate the signals between the 
science waveband and the waveband us ed for metrology; 

- a Mo dified Mach Zehnder (MMZ, ISerabvn & Colavital 
2001) beam combiner; 

- coupling optics to focus the outputs of the cross-combiner 
into single mode fibres. Chalcogenide fibre s can cover suc- 
cessfully the wavelength range 6.0-12 /zm dKsendzov et al.l 
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Fig. 2. Block diagram of the Darwin/TPF optical layout. Feed- 
back signals driving the tip-tilt/OPD control are represented by 
dashed lines. 



|2007|) and silver halide fibres can be used for modal filtering 
in at least the 10.5-17.5 /urn spectral range (further investiga- 
tions are however necessary to demonstrate that they are us- 
able i n the 17.5-20.0 /jm wavelength range, Ksendzov et al. 

120081) : 

- a detection assembly controlled at a temperature of 8 K and 
connected to the fibres. 

Darwin/TPF would be placed around the second Lagrange 
point (L2) by an Ariane 5 ECA vehicle. L2 is optimal to achieve 
passive cooling of the collector and beam combiner spacecraft 
down to 40 K by means of sunshades. An additional refrigerator 
within the beam combiner spacecraft cools the detector assem- 
bly to 8 K. Due to the configuration of the array and the need for 
solar avoidance, the instantaneous sky access is limited to an an- 
nulus with inner and outer hal f-angles o f 46° and 83° centred on 
the anti-sun vector (see Fi g . [Tl ICarld2005l) . This annulus transits 
the entire ecliptic circle during one year, giving access to almost 
the entire sky. 



2.2. Scientific objectives 

The main scientific objectives of Darwin/TPF are the detection 
of rocky planets similar to Earth and the spectroscopic analysis 
of their atmospheres at mid-infrared wavelengths (6 - 20 /mi). 
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Table 2. Parameters adopted for the performance simulations. 



Parameter 


Configuration 


Emma X- Array Design with a 6: 1 aspect ratio (see TableUJ 


Planet 


Earth-sized planet with a constant temperature of 265 K across the habitable zone 


Exozodiacal density 


Three times our solar system, based on Kelsall model (Kelsall et al. 1998) 


Mission duration 


5 years (2 years for survey and 3 years for characterisation) 


Retargeting time 


6 hours 


Integration efficiency 


70% (accounting for overhead loss) 


Target stars 


Darwin All Skv Survev Catalogue (DASCC, Kaltenegser et al. 2008) 


Earth-like planet per star (773,) 


1 


Habitable zone 


0.7- 1.5 AU scaled withLi /2 


Time allocation 


10% F, 50% G, 30% K and 10% M stars 


Array rotation period 


50000 seconds (not scaled with the array baseline length) 


SNR threshold for detection 


5 


SNR threshold for spectroscopy 


5 (C0 2 ,0 3 ) and 10 (C0 2 ,0 3 ,H 2 0) 



In addition to presenting an advantageous star/planet contrast, 
this wavelength range holds several spectral features relevant for 
the search of biological activity (CO2, H2O, and O3). The ob- 
serving scenario of Darwin/TPF consists of two phases, detec- 
tion and spectral characterisation, whose relative duration can 
be adjusted to optimise the scientific returrfj During the detec- 
tion phase of the mission (nominally 2 years), Darwin/TPF will 
examine nearby stars for evidence of terrestrial planets. A du- 
ration of 3 years is foreseen for the spectroscopy phase, for a 
total nominal mission lifetime of 5 years. An extension to 10 
years is possible and will depend on the results obtained during 
the 5 first years. Such an extension could be valuable to observe 
more M stars (only 10% of the baseline time being attributed to 
them), search for big planets around a significantly larger sam- 
ple of stars, and additional measurements on the most interesting 
targets already studied. 

The Darwin/TPF target star list has been generated from the 
Hipparcos catalogue, considering several criteria: the distance 
(< 30 pc), the brightness (< 12 V-mag), the spectral type (F, G, 
K, M main sequence stars), and the multiplicity (no companions 
within 1"). The corresponding star catalogue contains 1229 sin- 
gle main sequence st ars of which 107 are F, 23 5 are G, 536 are K, 
and 351 are M type (Kaltenegger et al. 2008). The survey of the 
Darwin/TPF stars and the possible detection of terrestrial plan- 
ets will start a new era of comparative planetology, especially by 
studying the relationship between habitability and stellar char- 
acteristics (e.g. spectral type, metallicity, age), planetary system 
characteristics (e.g. orbit), and atmospheric composition. 

3. Simulated performance 

3.1. The science simulator 

The performance predictions presented in this paper have been 
computed using the Darw in science SIMu lator developed at 
ESA/ESTEC (DarwinSIM. Iden Hartogll2005bl) . This simulator 
has been subject to extensive validation the past few years and 
its performance predictions were recently reconciled with a sim- 
ilar mission simulat or developed ind ependently at NASA/JPL 
for the TPF mission (lLav et alJ [2007). The two simulation tools 
have shown a very good agreement in SNR, giving similar in- 
tegration times for all the Darwin /TPF targets with a discrep- 
ancy lower than 10% in average dDefrere et al. 2008b). These 
two simulators have the same basic purpose. For a given instru- 
mental configuration and target catalogue, they assess the num- 

1 The detection phase might not be necessary if the targets are iden- 
tified in advance by radial velocity or astrometric surveys. 



ber of terrestrial planets that can be detected in the habitable 
zone of nearby main sequence stars and the number of possible 
follow-up spectroscopic observations during a nominal mission 
time. The duration of detection and spectroscopy phases can be 
adjusted to optimise the scientific return and is nominally set to 
2 and 3 years respectively. The parameters and assumptions used 
are summarized in Table [2] 

3.1.1. Detection phase 

The starting point of the simulations is the target star catalogue. 
Given a specific interferometer architecture, the simulator first 
identifies the stars which are observable from L2. For each of 
these observable stars, the basic calculation consists of an as- 
sessment of the required integration times to achieve an user- 
specified SNR for broad-band detection of a hypothetical Earth- 
like planet located inside the habitable zone. The habitable zone 
is assumed to be located between 0.7 and 1.5 AU for a G2V star 

1 II 

and is scaled with the square root of the stellar luminosity (L + ' ). 
Since the location of the planet around the star is a priori un- 
known, the integration time is computed from the requirement 
that it should ensure the detection of a planet for at least 90% 
of the possible locations in the habitable zone. Assuming plan- 
ets uniformly distributed along habitable orbits, this requires the 
computation of the probability distribution for finding a planet 
at a certain angular distance from the star. For each planetary 
position, the maximum SNR is computed by optimisation of the 
baseline length. The thermal flux of the habitable planet is as- 
sumed to be identical to that of Earth irrespective of the distance 
to the star. The exozodiacal clouds are simulated by assuming 
a nominal dust density 3 times larger than that in the solar sys- 
tem and a dust sublimation temperature of 1500 K. Under the as- 
sumption that the exozodiacal emission is symmetric around the 
target star, it will be suppressed by phase chopping, and therefore 
only contributes to shot noise. The noise sources included are the 
shot noise contributions from stellar leakage, local and exozo- 
diacal clouds, and instrumental infrared background. Instability 
noise is also present and is partly mitigated by phase chopping. 
A complete list can be found in Table [3] 

After the initial integration time assessment for detection, the 
targets are sorted by ascending integration time, removing from 
the list the targets for which the total integration time exceeds the 
total time during which they are visible from L2. Considering a 
slew time for re-targeting (nominally 6 hours) and an efficiency 
for the remaining observing time of 70%, the sorted list is cut 
off at the moment when the cumulative integration time exceeds 
the nominal survey period. Accounting for a specific time alio- 
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cation for each spectral type (10% F, 50% G, 30% K and 10% M 
stars), the resulting list defines the number of targets that can be 
surveyed during the detection phase. The actual number of plan- 
ets found will then depend on the number of terrestrial planets 
present in the habitable zone of target stars (77®). 

3.1.2. Spectroscopy phase 

The number of targets which can be characterised by spec- 
troscopy in a given time is computed similarly. The difference 
with the detection phase is that the integration times are com- 
puted for a given position in the habitable zone. The proper 
procedure would be to take into account all possible positions 
for the planet in a similar way to the detection phase but this 
would be far too time consuming. The strategy is then to con- 
sider only the most likely angular separation. Then, the total 
integration time is determined by the requirement to detect the 
absorption lines of O3, CO2 and H2O to a specified SNR. For 
the spectroscopy of CO2 and O3 (without H2O), an SNR of 
5 wo uld actually be sufficient for a secure detection (Fridlund 
I2005I) . Considering the spectroscopy of H2O is relatively more 
complex. Recent results suggest that, using a spectral resolution 
greater than 20, an SNR of 10 from 7.2 to 20 fim would be suf- 
ficient for H2O, CO2, and O3 spectroscopy (private communica- 
tion with F. Selsis, L. Kaltenegger and J. Paillet). In particular, 
these results suggest that the H2O band located below 7.2/im, 
which is much more time-consuming than the H2O band beyond 
17.2 fim, could be discarded. 

Two types of spectroscopic analysis are considered: the star- 
ing spectroscopy, where the array is kept in a position such that 
the planet resides on a peak of the modulation map, and the ro- 
tating spectroscopy, where the array keeps on rotating with re- 
spect to the target system so that the planet moves in and out of 
the peaks of the modulation map. Staring spectroscopy is more 
efficient in terms of signal acquisition, but requires an accurate 
knowledge of the planetary orbit. As for the detection phase, the 
total spectroscopy integration time should not exceed the total 
time during which the target is visible during the characterisa- 
tion phase. Accounting again for a given fraction of overhead 
loss, the targets are sorted with respect to ascending integration 
time, terminated where the cumulative time exceeds the length 
of the nominal characterisation period. The number of planets 
that can be characterised is then given with the assumption that 
there is one terrestrial planet in the habitable zone of each target 
star (77e=l). 

3.2. Coupling efficiency 

Coupling the optical beams into optical fibers is an essential 
part of the wavefront correction process, which is required for 
deep nulling. The theoretical efficiency of light injection into an 
optical fiber depends on several parameters: the core radius of 
the fiber, its numerical aperture, the wavelength, the diamete r 
of the telescope and its focal length (Ruilier & Cassaing 2 0011) . 
The method used in the simulator consists in choosing first the 
core radius of the fiber so as to ensure single-mode propaga- 
tion over the whole wavelength range. The f-number of the cou- 
pling optics can then be optimised to give the maximum cou- 
pling efficiency at a chosen wavelength and more importantly, to 
provide a roughly uniformly high coupling efficiency across the 
whole wavelength band. The coverage of the full science band 
of Darwin/TPF with one optical fiber is generally prohibited 
since the coupling efficiency drops rapidly with respect to the 
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Fig. 3. Coupling efficiency for Darwin/TPF with respect to the 
wavelength for an on-axis source and for a source with a fixed 
off-axis angle, corresponding to an Earth orbit around a Sun at 
10 pc. The core radius is chosen so as to stay single-mode on the 
whole wavelength. 



wavelength. Increasing the number of fibers improves the cou- 
pling efficiency but at the expense of complexity. In particular, 
it has been shown that the loss of targets for detection and spec- 
trosco py is about 5% bet ween the optimised 2-band and 3-band 
cases dden Hartogll20 05a). Considering the use of two fibers (re- 
spectively on the 6.0-12 /zm and 12-20 yum wavelength ranges), 
Fig.|3]shows the coupling efficiency for an on-axis source and for 
a source with a fixed off-axis angle, corresponding to an Earth or- 
bit around a Sun at lOpc. The coupling efficiency remains above 
70% over the whole wavelength band of each fiber. 

3.3. Modulation efficiency 

The modulation of the planetary signal during the observation is 
a direct consequence of the chopping process which is manda- 
tory to get rid of background noise sources such as exozodiacal 
and local zodiacal cloud emission. For the X-array configura- 
tion, the outputs of two Bracewell interferometers are combined 
with opposite phase shifts (± n/2) to produce two transmission 
maps (or "chop states"). Differencing the two transmission maps 
gives the chopped response of the interferometer, the so-called 
modulation map, which contains positive and negative values 
by construction (see Fig. |4}. Since the value of the modulation 
map varies across the field-of-view, the position and flux of the 
planet cannot be unambiguously inferred and an additional level 
of modulation is mandatory. This is provided by the rotation of 
the interferometer (typically with a period of 1 day). The plane- 
tary signal is therefore modulated as shown on the right part of 
Fig. [4] In order to retrieve the planetary signal, the most com- 
mon approach is correlation mapping, a technique closely re- 
lated to the Fo urier transform used for standard image synthe- 
sis CayJ[2005). The result is a correlation map, displayed for a 
single point source in the lower right part of Fig. [4] This repre- 
sents the Point Spread Function (PSF) of the array. This process, 
illustrated here for a single wavelength, is repeated across the 
wavelength range, and the maps are co-added to obtain the net 
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Fig. 4. Overview of phase chopping for the X-array configuration. Combining the beams with different phases produces two conju- 
gated transmission maps (or chop states), which are used to produce the chopped response. Array rotation then locates the planet by 
cross-correlation of the modulated chopped signal with a template function. 



correlation map. The broad range of wavelengths planned for 
Darwin/TPF greatly extends the spatial frequency coverage of 
the array, suppressing the side lobes of the PSF. 

After chopping and rotation, the part of the incoming sig- 
nal which is actually modulated and retrievable by synchronous 
demodulation is proportional to the "rotational modulation ef- 
ficiency". It is shown for the X-array configuration in Fig. [5] 
It depends on the radial distance from the star and reaches a 
peak value of 0.56 with an asymptotic value of 0.44. Since the 
planet position inside the habitable zone is a priori unknown, it 
is desirable that the effective modulation efficiency is as uniform 
as possible across the habitable zone to avoid too many recon- 
figurations of the interferometric array. Note that the rotational 
modulation effi ciency for several array configurations has been 
investigated bv lLavl d2005l) . 

3.4. Instability noise 

Instability noise is defined as the component of the demodulated 
nulled signal that arises from phase, amplitude and polarisation 
errors Ca3[2004). The power spectra of these instrumental ef- 
fects mix with each other so that perturbations at all frequencies, 
including DC, have an effect. Spacecraft vibrations, fringe track- 
ing offset, control noise, longitudinal chromatic dispersion, and 
birefringence are at the origin of the phase errors whereas tip/tilt, 
defocus, beam shear, and differential transmission produce am- 
plitude errors. These phase and amplitude errors induce a time- 
dependent asymmetry between the two chop states so that the 
modulation map does not remain centered on the nominal posi- 
tion of the line of sight (i.e. the position of the star). Hence a frac- 
tion of the starlight survives the modulation process and mixes 
with the planet photons. Although a simple binary phase chop 
removes a number of these systematic errors, it has no effect on 
the dominant amplitude-phase cross terms and on the co-phasing 
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Fig. 5. Rotational modulation efficiency for the Emma X-array 
with a 6:1 aspect ratio. 



errors. There is no phase chopping scheme that can remove the 
systematic errors without also removing the planetary signal. 

Three independ ent studies (Lay] 120041: Id'Arciol 120051: 
IChazelas et al.l 120061) have reviewed the instrumental require- 
ments on the Darwin/TPF mission that reduce the instrumental 
stellar leakage to a sufficiently low level for Earth-like planet 
detection. Assuming the presence of 1/f-type noise, these stud- 
ies showed that the requirements on amplitude and phase control 
are not driven by the null-floor leakage, but by i nstability noise. 
Considering a Dual-Chopped Bracewell (DCB, lLavll2004l) with 
4-m aperture telescopes operating at 10 yum, the different anal- 
yses show that a null depth of ~ 10~ 5 is generally sufficient to 
control the level of shot noise from the stellar leakage, but that a 
null depth of ~ 10~ 6 is required to prevent instability noise from 
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Table 3. Detailed view of the various contributors to the noise budget, given in photo-electrons persecond over the 6-20 yum wave- 
length range for Darwin/TPF in the Emma X-array configuration. The final SNRs are computed over a single rotation of 50000 s 
and for four targets representative of the target catalogue (see main text for further information). 
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Planetary signal [e-/s] 
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Photon noise [e-/s] 
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- Geometric leakage [e-/s] 
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0.019 
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0.008 


0.012 


- First order phase term [e-/s] 


0.021 


0.008 


0.002 


0.001 


- Second order phase-amplitude [e-/s] 


0.014 


0.012 


0.007 


0.012 


Total noise [e-/s] 


0.460 


0.388 


0.358 


0.365 


Integrated SNR (one rotation) 


16.0 


4.4 


2.1 


1.1 


Integration time for SNR=5 [h] 


1.36 


17.9 


81.8 


285 



becoming the dominant source of nois^J In particular, a 10~ 6 
null requires rms path control to within about 1.5 nm, and rms 
amplitude control of about 0.1%. 

In order to relax these very str ingent requirements, several 
techniques have b een investigated (Lay 20051 lLane et al.ll2006t 
Gabor et al. 2008). Discussion of these mitigation techniques is 
beyond the scope of this paper, where we assume that instability 
noise is sufficiently low to ensure the H2O spectroscopy at 7-yizm 
of an Earth-like planet orbiting aro und a Sun located at 15 pc. 
Applying the analytical method of (Lav 2004) to the Emma X- 
array with the parameters listed in Table [TJ we derive the con- 
straints on the instrument stability such that instability noise is 
dominated by shot noise by a factor 5 at 7 fim over one rota- 
tion period (f lot ) of 50000 seconds (with a spectral resolution 
of 20). Considering 1/f-type PSDs defined on the [l/f lot ,10 4 ] Hz 
range, this corresponds to residuals rms OPD and amplitude er- 
rors of about 1.5 nm and 0.05% respectively. These values will 
be used thorough this study (see appendixlAlfor further details). 
Although our computation has been done at 7 //m where insta- 
bility noise is much higher than at 10 /zm, the se constraints are 
not far from the values derived by lLavl d2004t) for two reasons. 
First, the telescopes considered here are smaller so that shot 
noise is relatively more dominant than in the previous analyses 
(shot noise is proportional to the square root of the stellar flux 
while the planetary signal and instability noise are directly pro- 
portional to the stellar flux). Secondly, the interferometer con- 
figuration is stretched by comparison with the DCB so that the 
planetary signal is modulated at higher frequencies, where the 
instability noise is lower assuming 1/f-type PSDs. 

3.5. Signal-to-noise analysis 

In this section, we present the different sources of noise simu- 
lated by DarwinSIM and the level at which they contribute to the 
final SNR for four targets representative of the Darwin/TPF cat- 

2 The current state-of-the-art for broadband nulling experiments 
is a 10~ 5 null which has been recentl y demonstrated a t 10 fim 
(A/I/lambda=34%) with the adaptive nuller dPeters. et ai1l2009h . 



alogue: an MOV star located at 5 pc, a K0V star at lOpc, a G2V 
star at 15 pc and an F0V star at 20 pc. The noise budget of each 
source is shown in Table [3] for a single rotation of 50000 s and 
for the optimum baseline length (computed by minimizing the 
integration time). The different contributors are described here- 
after. 

- The stellar signal represents the total number of photo- 
electrons that are generated by stellar photons detected in 
both constructive and destructive outputs. 

- The planetary signal is the demodulated amount of photo- 
electrons that come from an Earth-like planet located at 1 AU 
from the star. 

- Shot noise is due to the statistical arrival process of the 
photons from all sources. It comes from the contributions 
from stellar leakage, the exozodiacal dust, the local zodia- 
cal cloud, the thermal emission from the telescopes, and the 
stray light. 

- Geometric stellar leakage accounts for the imperfect rejec- 
tion of the stellar photons due to the finite size of the star 
and the non-null response of the interferometer for small off- 
axis angles. 

- Null-floor leakage accounts for the stellar photons that leak 
through the output of the interferometer due to the influ- 
ence of instrumental imperfections such as co-phasing er- 
rors, wavefront errors or mismatches in the intensities of the 
beams. 

- The 3-zodi signal is the shot noise contribution from the 
circumstellar disc, assumed to be f ace-on and to follo w the 
same model as in the solar system dKelsall et al Jl 19981) . ex- 
cept for a global density factor of 3. 

- Local zodiacal signal is the shot noise contribution from the 
solar zodiacal cloud, taking into account the spacecraft loca- 
tion at L2 and the pointing direction. 

- Thermal background accounts for the emission of the tele- 
scopes. 

- Stray light is made of the photons originating from outside 
the interferometer and which do not follow the nominal route 
to the detector. It includes scattered light from the target star, 
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Fig. 6. Input (a) and detected (b) signals for an Earth-like planet 
orbiting at 1 AU around a G2V star located at 15 pc. The demod- 
ulated signals are computed over a single rotation of 50000 s. 



thermal photons from the instrument and any solar photon 
that are scattered into the instrument. We assume a nominal 
value of 10 photons per second and per spectral channel. 

- Dark current is the constant response produced by the detec- 
tor when it is not actively being exposed to light. We con- 
sider a nominal value of of 4 electrons rms per read and per 
spectral channel. 

- Detector noise is computed assuming a read-out noise of 4 
electrons rms and a typical read-out frequency of 1 Hz. 

- Instability noise has been discussed in section [3741 It is com- 
puted for rms OPD and amplitude errors of 1.5 nm and 
0.05% respectively (and defined on 1/f-type power spectra). 

As expected, shot noise is the dominant contributor for all 
targets (we have derived the constraints on the instrument sta- 
bility so that instability noise is not dominant in section 13.41 ). 
The emission of the local zodiacal cloud is generally the dom- 
inant source of shot noise but geometric stellar leakage and the 
exozodiacal dust also produce a significant part of it. Instability 
noise is dominated by s econd order phase-amplitude cross terms 
as pointed out by lLavl {2004) and contributes weakly to the in- 
tegrated SNR. The dependence on the wavelength is shown in 
Fig. [6] for the G2V star located at 15 pc. The upper figure rep- 
resents the different input signals, showing that stellar leakage 
dominates at short wavelengths while the local zodiacal cloud 



Table 4. Expected performance in terms of number of stars sur- 
veyed and planets characterised during the nominal 5-year mis- 
sion for various telescope diameters and planet radii. All stars 
are assumed to host a planet in the habitable zone and to be sur- 
rounded by an exozodiacal cloud 3 times denser that in the solar 
system. 
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emission is dominant at long wavelengths. The lower figure 
represents shot noise, instability noise and the detected signal 
from the planet. The SNR is also represented and measured on 
the right-hand side vertical axis. It is maximum around 10/im, 
where the planetary signal peaks, and decreases rapidly toward 
short wavelengths, suggesting that the spectroscopy at these 
wavelengths would be the most time-consuming. Integrating the 
signal and noise sources from 6 to 20 fim gives an SNR of about 
2 (see table |3]l so that about 6 rotations of the interferometric 
array would be necessary to achieve the detection of the planet 
with an SNR of 5. For spectroscopy, the integration time has to 
be much longer since the individual SNR in each spectral chan- 
nel is significantly lower. For instance, about 150 rotations of the 
interferometric array would be necessary to achieve an SNR of 
10 at 10/im. 

3.6. Expected performance 

Considering the assumptions given in Table 12 the simulated 
performance of Darwin/TPF is shown in Table [4] for various 
aperture sizes and planet radii. Considering Earth-radius plan- 
ets within the habitable zone, about 200 stars, well spread among 
the four selected spectral types, can be surveyed during the nom- 
inal 2-year detection phase. This number reaches about 500 with 
4-m aperture telescopes. Darwin/TPF will thus provide statisti- 
cally meaningful results on nearby planetary systems. As already 
indicated by Table[3] nearby K and M dwarfs are the best-suited 
targets in terms of Earth-like planet detection capabilities. 

For the spectroscopy phase, a required SNR of 5 has been 
assumed for the detection of CO2 and O3, as discussed in sec- 
tion |3.1| For the full characterisation (i.e. searching for the pres- 
ence of H2O, CO2, and O3), the required SNR has been fixed to 
10 on the 7.2-20-^m wavelength range. With these assumptions, 
CO2 and 03 could be searched for about 40 planets (resp. 60) 
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Fig. 7. Impact of the exozodiacal dust density on the SNR for 
different target star s. The exozodiacal disc is assumed to follow 
the Kelsall model (Kels all et alj[l9 98) and to be seen in face-on 
orientation. The horizontal dotted line corresponds to an increase 
of integration time by a factor 2 with respect to the 0-zodi case. 
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Fig. 8. Maximum number of zodis with respect to the target dis- 
tance for the Darwin/TPF target stars. The maximum number of 
zodis corresponds to an increase of integration time by a factor 
two with respect to the 0-zodi case. 



with rotational spectroscopy (resp. staring spectroscopy) while 
H2O could potentially be detected on 20 (resp. 30) planets 
during the 3-year characterisation phase. These values would 
be roughly halved for 1-m aperture telescopes. Although star- 
ing spectroscopy presents (as expected) better results, rotational 
spectroscopy is more secure since it does not rely on an accu- 
rate localisation of the planet. It is also interesting to note that in 
the case of planets with radii 1.5 time as large as that of Earth, 
the number of planets for which H2O spectroscopy could be per- 
formed is doubled. 

4. Impact of the exozodiacal cloud density 

The amount of exozodiacal dust in the habitable zone of nearby 
main sequence stars is one of the main design drivers for 
the Darwin/TPF mission. Depending on their morphology and 
brightness, exozodiacal dust clouds can seriously hamper the 
capability of a nulling interferometer to detect and characterise 
habitable terrestrial planets. Under the assumption that it is cen- 
trally symmetric around the target star, the exozodiacal cloud 
is suppressed by the chopping process, and therefore only con- 
tributes to shot noise. An exozodiacal cloud similar to the local 
zodiacal disc emits 350 times as much flux at 10/im than an 
Earth-like planet, so that it generally drives the integration time 
as the disc becomes a few times denser than the local zodiacal 
cloud. A previous study performed for the DCB with 3-m aper- 
ture telescopes observing a G2V star located at lOpc has led 
to the conclusion that detecting Earth-like planets around a star 
for which the exoz odiacal cloud density is larger than 20 zodis 
would be difficult dBeichman et alj|2006l) . Nevertheless, the tol- 
erable amount of dust around a nearby main sequence star highly 
depends on several parameters such as the telescope size, the 
target distance and spectral type. Another parameter which can 
affect the performance of the interferometer is the presence of 
asymmetric structures in the exozodiacal disc such as clumps or 
offset due to the presence of planets. These asymmetries have a 
different impact on the mission performance because they intro- 
duce a signal which is not perfectly suppressed by phase chop- 
ping and can mimic the planetary signal. This section is focused 
on the impact of the exozodiacal dust density on the integration 



time and the consequence on the number of targets that can be 
surveyed during the mission lifetime. The impact of asymmetric 
structures is discussed in section|5] 

4. 1 . Analysis per individual target 

Currently, very little is known about the amount of exozo- 
diacal dust in the habitable zone of nearby main sequence 
stars. First results have been obtained only very recently using 
classical infrared interferometry at the CHARA array (Mount 
Wilson, USA) and at the VLTI (Cerro Paranal, Chile). These 
instruments have revealed the presence of hot dust in the in- 
ner part of planetary systems around a few nearby main se- 
quenc e stars with a sensitivity of approximately one thousand 
zodis dAbsil et alj2006tlDi Folco et alj2007tlAbsil et"aill2008bl 
2009: lAkeson et alj|2009l) . Recent observations using ground- 
based nulling interferometry at the Keck observatory (Mauna 
Kea, USA) have shown impro ved sensitivity to e xozodiacal dust 
clouds of a few hundred zodis dStark et alj |2009). Given the lack 
of information on exozodiacal clouds with densities of a few 
zodis, we investigate in this section the impact of exozodiacal 
dust density on the performance of Darwin/TPF. 

Considering centrally symmetric face-on exozodiacal discs, 
Fig. Q shows the SNR (normalised to the SNR for no exozodi- 
acal cloud) integrated over the 6-20 fim wavelength range as a 
function of the exozodiacal dust density for the 4 typical target 
stars used in section I3.5I We consider the normalised SNR be- 
cause it does not depend on the integration time (the planetary 
signal is removed from the equation) which has the advantage 
to provide a common basis to compare the different target stars. 
Looking at Fig. [7] the impact of the exozodiacal dust density on 
the normalised SNR is particularly harmful for the hottest tar- 
get stars which present the brightest exozodiacal discs. For the 
F0V star located at 20 pc, the normalised SNR is reduced by 
a factor of about 2.5 between the and 100-zodi cases while 
for the MOV star located at 5 pc it is only reduced by a factor 
of about 1.5. Assuming that the integration time should not be 
twice longer than in the 0-zodi case (see the horizontal dotted 
curve), the maximum number of zodis are about 70, 50, 50 and 
20 respectively for the MOV star located at 5 pc, the K0V star at 
lOpc, the G2V star at 15 pc and the F0V star at 20 pc. 
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Fig. 9. Tolerable exozodiacal dust density as a function of the 
number of targets that can be observed during the mission life- 
time normalised to the zodi-free case (R bsX accounting for over- 
heads and a specific time allocation for each spectral type (10% 
F, 50% G, 30% K and 10% M). 



The distance of the star also plays an important role. As the 
distance to the target system increases, the flux collected from 
the exozodiacal cloud decreases while the flux collected from 
the local zodiacal cloud remains the same for all targets. The 
contribution of the exozodiacal dust cloud to the noise level be- 
comes therefore relatively less important so that a higher dust 
density can be tolerated around the target. This explains why the 
curves of the G2V star and the K0V star almost coincide despite 
the fact that the G2V star is hotter. This behavior is illustrated 
in Fig. [8] showing the maximum number of zodis with respect 
to the target distance for the whole Darwin/TPF catalogue. This 
maximum number of zodis corresponds to an increase of inte- 
gration time by a factor two with respect to the 0-zodi case. It 
depends on the target distance and spectral type, and can take a 
value from few zodis up to several hundred zodis for the most 
distant stars. For a given distance to the target system, the maxi- 
mum number of zodis increases with the stellar temperature, the 
zodi constraint being more severe on F stars than on M stars, 
while for a given spectral type, the zodi tolerance increases with 
the target distance. 

4.2. Tolerable dust density 

So far we have examined the maximum exozodiacal dust den- 
sity for each target which corresponds to an integration time per 
target equal to twice the zero-zodi integration time. However, 
this does not tell us anything about the total number of stars that 
could be observed over the survey time of the mission. Here we 
derive the exozodiacal dust density that can be tolerated around 
nearby main sequence stars so that a given number of stars can 
be observed during the nominal mission lifetime. 

To calculate the total number of targets that can be surveyed 
during the mission lifetime, we compute the integration time for 
each observable target as a function of the exozodiacal dust den- 
sity and add them in ascending order as described in section [37X1 
Considering a slew time of 6 hours and an efficiency for the 
remaining observing time of 70%, the list is cut off when the 
cumulative integration time exceeds the nominal survey period. 
Applying this procedure to each spectral type (with the corre- 
sponding time allocation, see Table [2]) and to the whole target 



list, Fig.|9]shows the tolerable dust density as a function of R b s , 
the ratio of the number of targets that can be observed during 
the mission lifetime in the presence of exozodiacal clouds of the 
given density to the number of targets that can be observed dur- 
ing the mission lifetime in the absence of exozodiacal clouds. 
The corresponding number of target stars that can be observed 
can easily be computed using Table [4] 

As exozodiacal discs become denser, R Q b s decreases so that 
the tolerable dust density depends on the goal of the mission in 
terms of the number of stars that have to be observed. In order 
to survey approximately 50% of the stars that could be observed 
in the absence of exozodiacal discs, a dust density as high as 
400 zodis can be tolerated, while only 30 zodis are tolerable to 
observe 90% of the stars. Considering that at least 150 target^ 
have to be observed during the mission lifetime (about 75% of 
the stars), exozodiacal discs with a density of about 100 zodis 
can be tolerated around the targets stars. This tolerable dust den- 
sity is computed for the whole catalogue assuming that the time 
allocation on each spectral type is maintained. In practice, the 
effect of exozodiacal dust is more pronounced for early type 
stars. This behavior is illustrated in Fig.|9]for each spectral type. 
For a dust density of 100 zodis, about 40%, 60%, 80% and 90% 
of the F, G, K and M stars respectively can still be surveyed. 
Conversely, in order to survey at least 75% of the stars, the tol- 
erable dust densities are about 10, 50, 100 and 300 zodis for F, 
G, K and M stars respectively. 

These results show how important it is to observe in ad- 
vance the Darwin/TPF targets in order to maximize the num- 
ber of stars that can be surveyed during the mission lifetime. 
Ground-based nulling instru ments like LB TI ("Large Binocular 
Telescope Interferometer", iHinz et al.l l2008h and ALADDIN 
('Antarctic L-band Astr ophysics Discovery Demonstrator for 
Interferometric Nulling". lAbsil et aDl2008ah would be ideal to 
reach the detection of 50-zodi exozodiacal discs with a sky cov- 
erage sufficient to observe almost the entire Darwin/TPF cata- 
logue. 

4.3. Influence of the telescope size 

Increasing the telescope diameter has different influences on the 
individual signal and noise sources. The planetary signal in- 
creases as D 2 while the shot noise contributions from geometric 
leakage and exozodiacal cloud increase as D. The relative con- 
tribution from the local zodiacal cloud to shot noise is reduced 
for larger aperture telescopes, due to the smaller field-of-view. 
Since the local zodiacal cloud emission is generally one of the 
dominant noise sources, the relative impact of the exozodiacal 
cloud density on the SNR becomes therefore more significant 
for larger telescopes. Using the assumptions of Tabled this be- 
havior is illustrated in Fig. [10] which shows the tolerable exo- 
zodiacal dust density with respect to R Q b s for different aperture 
sizes. 

For a given dust density, the loss of observable targets dur- 
ing the mission lifetime with respect to the case without exo- 
zodiacal disc is more important for larger aperture telescopes. 
For instance, the tolerable densities are 50, 30 and 15 zodis re- 
spectively for 1-m, 2-m and 4-m apertures telescopes in order to 
survey at least 90% of the nominal targets. These values become 
200, 100 and 60 zodis if 75% of the targets have to be surveyed. 
Considering again that 150 targets have to be observed during 
the mission lifetime, dust densities as high as 100 and 600 zodis 



3 Detecting 150 targets has been defined by both ESA and NASA as 
the minimum mission requirement for Darwin/TPF-1. 
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Fig. 10. Tolerable exozodiacal dust density for different aper- 
ture sizes as a function of the number of targets that can be ob- 
served during the mission lifetime normalised to the zodi-free 
case (Robs)- The results are computed for the nominal mission 
lifetime accounting for overheads and a specific time allocation 
for each spectral type (10% F, 50% G, 30% K and 10% M). 



can be tolerated around the target stars for 2-m and 4-m aper- 
tures telescopes respectively (150 targets being not detectable 
within the survey time with 1-m aperture telescopes). Due to 
the better nominal performance achieved with 4-m aperture tele- 
scopes (about 500 targets surveyed during the survey time, see 
Table HI, the maximum dust density to survey 150 targets (600 
zodis) is higher than for 2-m aperture telescopes (100 zodis) but 
the corresponding loss of surveyed targets with respect to the 
nominal case is much more important (R ODS of 30% vs 75%). 
In practice, larger aperture telescopes would obviously be better 
to maximize the scientific performance but the final choice will 
also result from a trade-off with cost and feasibility. 

5. Impact of the exozodiacal cloud morphology 

The previous results have been obtained assuming that the 
exozodiacal dust emission is centrally symmetric around the 
target star so that it is suppressed by phase chopping (and 
therefore only contributes to shot noise). However, exozodi- 
acal discs are likely to show resonant structures or an off- 
set with respect to the central star due to the the gravita- 
tional influence of embedded planets. Thes e resonant structures 
have been predicted by theoretical stud ies (Roques et al. 1994; 
iLiou & Zookll 19991: lOzernoy et al.ll206"ol) . and similar structures 
have been observed in few cases around nearby main-sequence 
stars (e.g.. IWilner et alj|2002t iGreaves et all 120051: l kalas et alj 
l2005t ISchneider et alJl2009l) . The most well-studied example of 
an asymmetric disc is the solar zodiacal cloud, which exhibits 
seve ral structures interpreted as t he dynamical sig nature of plan- 
ets (iDermott et all 1 19851 1 19941: iReach et al.llT995h . This trend 
suggests that exozodiacal clouds may be full of rings, clumps, 
and other asymmetries induced by the presence of embedded 
planets. These asymmetric structures around the target star are 
not perfectly canceled by the phase chopping process and part 
of the exozodi signal can then mimic the planetary signal. If 
the demodulated contribution from the exozodiacal disc is sig- 
nificantly higher than that of the planet, it would be difficult to 
isolate the pla netary sign al, whatever the integration time. As 
mentioned by lLayl d2004h . asymmetric inhomogeneities at the 



0.1% level of the total exozodiacal flux can be confused with 
a planetary signal. The problem becomes even more serious as 
the dust density increases: a 10-zodi exozodiacal disc must be 
smooth at the 0.01% level, in the region to be searched for plan- 
ets. Considering asymmetric exozodiacal discs, we derive in this 
section the tolerable dust density in order to ensure the detection 
of Earth-like planets with Darwin/TPF. 

5.1. Methodology 

At the output of the interferometer, the to tal detecte d photon rate 
(excluding stray light) can be written as (lLavll2004h : 

N= \ \ B sky (s)R(s)P(s)Gd8da, (1) 
Jg J a 

where B^y is the brightness distribution on the sky for a spec- 
tral channel centered on wavelength A and a bandwidth AA <<c A 
(units in photons/s/m 2 ), s is a unit vector whose direction repre- 
sents position on the sky, P(s) is a field-of-view taper function 
resulting from the size of a collecting aperture and from the re- 
sponse of the single-mode spatial filter, and R(s) is the intensity 
response of the interferometer on the sky (excluding the taper), 
the so-called transmission map. Expression Q] can be written as 
a function of the double Fourier transform of the sky brightness, 
including the field-of-view taper, which are denoted by a hori- 
zontal bar: 

N = ^ 2j [ COS (0; - <f>k) B S ky,jk,.sym 

j k 

- Sin((f>j - 0/t)fiskyjk,asym ], (2) 

where B s k y , sy m (resp. B s k y , a sym) is the symmetric (resp. asymmet- 
ric) part of the sky brightness distribution and <pj the phase re- 
sponse of the telescope j. The total photon rate is therefore a sum 
over all possible pairs of collectors. The baselines with a phase 
difference that is an integer multiple of n (0, +n, +2n) couple en- 
tirely to the symmetric brightness distribution (star, local zodia- 
cal cloud, exozodiacal disc). Any asymmetric brightness distri- 
bution couples entirely to baselines with a phase difference that 
is an odd multiple of n/2 (e.g., planet, clump). From Eq.|2j the 
detected photon rate from the exozodiacal disc is simply given 
by: 

NeZ = Zj AjAk [ C0S ^j ~ ^) B EZ,jk, S ym 
j k 

- Sin(0j - t )B E Z,jk,asym] , (3) 

The demodulated signal from the exozodiacal disc can then be 
obtained (after combination of the two chop states): 




where T is the integration time, N^ z (resp. N^ z ) the contribu- 
tion of the exozodiacal disc at the output of the left (resp. right) 
chop state and 77 the demodulation template function. This tem- 
plate functi on is used to extract t he planetary signal by cross- 
correlation (Angel & WooH)ll997l) . It represents the time series 
(normalized to have a rms value of one) that would be obtained, 
should a planet be present at a given position. Taking into ac- 
count that the phase response of the interferometer is such that 
<fi = -(pf with the 0f equal to (n/2, 0, n, 3tt/2), Eq.[3]and|4]give: 
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Fig. 11. Upper: Thermal flux (6-20 pm) produced by a 10-zodi exozodiacal dust cloud a round a G2V sta r for four different disc 
inclinations (0°, 30°, 60° and 90°) and assuming a Dohnanyi distribution of particle sizes (|Pohnanvill l969). The images have bee n 
simulated assuming a Earth-mass planet located at 1 AU on the x-axis (90 degrees clockwise from vertical. I Stark & Kuchnerl20 08). 
Lower: Corresponding asymmetric brightness distributions obtained by subtracting to each pixel its centrally symmetric counterpart. 



Cez = — V V A jA k sin(0 ; - - <f> k ) \ B EZ j k . RS y m T]dt , (5) 

which shows that only asymmetric components of the brightness 
distribution contribute to the demodulated signal. More specif- 
ically, only the baselines with a "fractional-^" phase difference 
can produce a demodulated signal. Such baselines are mandatory 
in order to generate odd harmonics of the rotation frequency. A 
solution to mitigate the influence of the asymmetric structures in 
the disc is to have long imaging ("fractional-^-") baselines that 
resolve out the more spatially extended emission from the exo- 
zodi variations, leaving only the point like emission from plan- 
ets. 

5.2. Impact of clumps 

The origin of some asymmetric clumpy structures, i.e. local 
density enhancements, in exozodiacal discs may be attributed 
to to the gravitational influence of planets on the small dust 
grains. After their release from parent bodies via collisions or 
outgassing, dust grains experience different paths in the stellar 
system, depending on their effective size. Whereas the small- 
est particles are ejected from the planetary systems by radiation 
pressure in a dynamical time, larger p articles slowly sp iral in- 
ward due to Pointing-Robertson drag (Rob ertsonlll937l) . While 
spiraling toward their host star, dust particles may become tem- 
porarily trapped in mean motion resonance with planets, extend- 
ing their lifetimes. This trapping locally enhances the particle 
density, creating structures, originally described for the solar 
zodiacal cloud as circumstellar rings, bands, and clumps (e.g., 
iKelsall et alJfl998h . 

In order to address the impact of such structures on 
the performance of D arwin/TPF, we use the results of 
Stark & Kuchnerl d2008l) . who synthesized images of cir- 
cumstellar discs with resonant rings structures due to em- 
bedded terrestrial-mass planets. Among the studies that 



have examined the geometry of these resonant signatures 
(e.gjKuchner & Holmaij2003ilMoro-Martfn & Malhotral2005l: 
Reche et al. 2008; Stark & Kuchner 2008), these images are par- 
ticularly convenient for our study since they include enough par- 
ticles to overcome the limitations of previous simulations, which 
were often dominated by various sources of Poisson noise, and 
allow for quantitative study of the modeled ring structures. In 
addition, these images are geared toward terrestrial-mass plan- 
ets at a few AU from the star, whereas most other studies con- 
cern more massive planets located much farther from the star. 
We used the lStark & Kuchnerl (120081) models to produce thermal 
emission images of inclined discs (0°, 30°, 60° and 90°) with 
resonant ring structures. We investigated disc models for a sys- 
tem with an Earth-mass planet on a circular orbit at 1 AU around 
a G2V star located at 15 pc and for a Dohnanyi distribution rang - 
ing in size from the blowout size up to 120 pm dDohnanvill 969). 
The thermal emission produced by such exozodiacal discs are 
given in the upper part of Fig. QT|in a wavelength range of 6- 
20 pm. 

The images given in Fig. QT| can be thought of as upper- 
limits to the brightness o f structures due to an Earth-like planet. 
IStark & Kuchnerl Q2008) ignored dust from parent bodies with 
large inclinations and eccentricities, such as comets, which 
would tend to wash out any resonant structure. Additionally, 
these models ignore the effects of collisions, which smooth out 
overdense regions of the disc and reduce azimuthal asymmetries 
(IStark & Kuchnerl l2009). In every simulation, the parent bodies 
were initially distributed from 3.5 to 4.5 AU in an asteroid belt- 
like ring. The Earth-mass planet is oriented along the x-axis (lo- 
cated at 66 mas on the x-axis) with a noticeable gap in the ring 
at its position. The models are truncated at half the semi-major 
axis of the planet, resulting in the inner holes in the images of 
Fig. Q~T] In reality, the dust density distribution should continue 
inward to the dust sublimation radius in the absence of additional 
perturbers. This "missing" inner disc should however not affect 
our results, since the inner disc would be centrally symmetric so 
that it does not contribute to the detected signal. 
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Fig. 12. Upper: Chopped photon rate from an Earth-like planet and a 10-zodi asymmetric disc with respect to the rotation angle for 
two disc inclinations (0° and 60°) and different wavelength ranges (at 10pm and in the full wavelength range). Lower: Corresponding 
dirty map formed from the cross correlation of the measured signal with templates of the signal expected from a point source at 
each location on the sky. 



Introducing these images into DarwinSIM, we compute the 
chopped photon rate from the exozodiacal disc as a function of 
the array rotation angle. This is represented in the upper part of 
Fig.[12]for two different disc inclinations (0° and 60°) at 10 pm 
(left figures) and for broadband detection (6-20/im, right fig- 
ures). The density of the disc has been scaled up to 10 zodis 
and the chopped planetary signal represented for comparison 
(dashed curve). In all cases, it is dominated by the chopped 
signal from the exozodiacal disc (solid curve) and particularly 
for high disc inclinations. To disentangle the planetary signal 
from the disc signal, it is necessary to apply the cross-correlation 
method to build the so-called dirty map. Applying cross corre- 
lation of the measured signal (disc + planet) with templates of 
the signal expected from a point source at each location on the 
sky (computed using Eq. [5]), the dirty maps represented in the 
lower part of Fig. [12] are obtained. This process transforms the 
rotationally modulated signal into a map of the sky by cross- 
correlation, which is equivalent to the Fourier transform used in 
standard synthesis imaging. The planet is located at 66 mas on 
the x-axis and presents a demodulated signal of about 0.028 e~/s 
at 10/im and about 0.62 e~/s for broadband detection. The to- 
tal demodulated signal at the planet position is however lower 
due to the negative contribution from the exozodiacal disc which 
presents a hole near the planet. The demodulation of the exozodi- 
acal discs also produces main peaks which are maximum around 
30-50 mas from the host star (in agreement with the asymmetric 
brightness distribution in the initial images, see the lower part of 
Fig.rTTTi. Fortunately, the high angular resolution provided by the 
long imaging baseline is sufficient to spatially distinguish these 
components from the planetary signal and only the contribution 
from the hole around the planet significantly contributes to the 
noise level. 

In order to ensure the planet detecti on, we adopt a crite- 
rion commonly used in AO imag ing dMacintosh et alJ 12003b 
iHinklev et aLll2007t ISerabvnll2009l) . The noise level is taken to 
be the rms deviation of the pixel counts within an annulus of 
width equal to the size of the PSF at half maximum. Considering 
a detection threshold of 5, the results are given in Table [5] for 
different wavelengths and disc inclinations. The tolerable disc 



density ranges between about 1 and 15 zodis, depending on the 
disc inclination. The detection is particularly difficult for highly 
inclined discs for which the asymmetric components are more 
dominant and at long wavelengths where the planetary signal is 
weaker. Combining the spectral channels to obtain a broadband 
correlation map (see Fig.fT2l) reduces the impact of sidelobes as- 
sociated with each main peak but does not significantly improve 
the results. This is because the performance are limited by the 
main peak induced by the hole near the planet rather than by 
sidelobes. Note that this hole generally induces a response 2 to 
4 times larger the rms deviation of the pixel counts within the 
annulus so that it might be marginally interpreted as a planet de- 
tection (a "false positive"). However, it could also be seen as 
an indirect way to detect the presence of a planet within the 
hole since such compact structure is expected to be created by 
a planet. 

In order to relax these stringent constraints on the exozo- 
diacal dust density, more sop histicated techniques of dat a pro- 
cessing could be useful (e.g jThiebaut & Mugnierl 12006). The 
capability of these techniques still needs to be investigated and 
is beyond the scope of this paper. The only secure way to en- 
sure the detection and characterization of Earth-like planets with 
DARWiN/TPF-like missions is to observe in advance the nearby 
main-sequence stars in order to remove from the target list the 
stars with a too high inclined/bright exozodiacal disc. To achieve 
this goal, a space-based nulling interferometer such as FKSI 
("Fourier-Kelvin Stellar Interferometer") would be ideal with a 
sensitivity sufficient to detect exozodiacal discs down to 1 zodi 
dDefrere et al.ll2008al) . 

5.3. Impact of the disc offset 

An offset between the center of symmetry of a dust cloud and 
its host star is a natural consequence of the gravitational interac- 
tion with planets. In the solar system, the center of the zodiacal 
cloud is shifted by about 0.013 AU from the Sun due mostly to 
Jupiter dLandgraf & Jehnll2001l) . The offset can be much larger, 
as sho wn in the case of t he Fomalhaut system with an offset of 
15 AU dKalas et al.ll2005l) . Even when inhomogeneities such as 
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Table 5. Tolerable exozodiacal dust density for different disc in- 
clinations and wavelengths. The detection threshold is taken to 
be 5 times the rms deviation of the pixel counts within an annu- 
lus of width equal to the size of the PSF at half maximum. 



Disc incl. 


8yum 


lOyum 


16yum 


Wide 


0° 


12.2 


15.3 


7.0 


14.0 


30° 


9.1 


15.0 


6.9 


13.8 


60° 


6.1 


6.8 


2.1 


3.8 


90° 


1.1 


1.4 


1.0 


2.6 



clumps are not present, an offset cloud produces an asymmetric 
brightness distribution such that a part of the exozodiacal disc 
signal survives the chopping process. Using the zodipic pack- 
ag^ we produce images of solar-like zodiacal discs with a given 
offset and use them to compute the demodulated signal at the 
output of the interferometer. The results are presented in Fig.QjJ 
showing the tolerable dust density with respect to the disc off- 
set for a G2V star located at 15 pc and for different wavelengths. 
The disc is assumed to be seen face-on. 

As the distance between the host star and the center of sym- 
metry of the exozodiacal disc increases, the tolerable dust den- 
sity to detect an Earth-like planet located at 1 AU becomes more 
severe and reaches less than 5 zodis for an offset of 1 AU. For in- 
dividual spectral channels, the tolerable dust density rapidly de- 
creases to reach the value of 20 zodis at 0.05, 0. 15 and 0.25 AU 
respectively at 8/im, 10 yum and 16 yum. The results for broad- 
band detection are much better with a tolerable dust density of 
20 zodis only for a disc presenting an offset larger than 0.6 AU. 
Considering a tolerable exozodiacal dust density of 100 zodis, 
the offset between the host star and the center of symmetry of 
the exozodiacal disc can be as high as 0.4 AU in order to ensure 
the planet detection. For an offset similar to that of the solar zo- 
diacal cloud (about 0.013 AU), this tolerable dust density is even 
much higher (few thousand zodis). 

6. Conclusions 

Infrared nulling interferometry is the core technique of future 
life-finding space missions such as ESA's Darwin and NASA's 
Terrestrial Planet Finder (TPF). Observing in the infrared (6-20 
/um), these missions will be able to characterise the atmosphere 
of habitable extrasolar planets orbiting around nearby main se- 
quence stars. This ability to study distant planets strongly de- 
pends on exozodiacal clouds around the stars, which can hamper 
the planet detection. Considering the nominal mission architec- 
ture with 2-m aperture telescopes, we show that centrally sym- 
metric exozodiacal dust discs about 100 times denser than the 
solar zodiacal cloud can be tolerated in order to survey at least 
150 targets during the mission lifetime. The actual number of 
planet detections will then depend on the number of terrestrial 
planets in the habitable zone of target systems. 

The presence of asymmetric structures in exozodiacal discs 
(e.g., clumps or offset) may be more problematic. While the 
cloud brightness drives the integration time necessary to dis- 
entangle the planetary photons from the background noise, the 
emission from inhomogeneities are not perfectly subtracted by 
phase chopping so that a part of the disc signal can mimic the 
planet. To address this issue, we consider modeled resonant 
structures produced by an Earth-like planet and introduce the 
corresponding image into DarwinSIM, the mission science sim- 
ulator. Even for exozodiacal discs a few times brighter than the 

4 http;//asd.gsfc. nasa.gov/Marc.Kuchner/liome. html 
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Fig. 13. Tolerable exozodiacal dust density with respect to the 
offset between the center of symmetry of the exozodiacal disc 
and the central star (a G2V star located at 15 pc). The disc is 
assumed to be seen in face-on orientation. 



solar zodiacal cloud, the contribution of these asymmetric struc- 
tures can be much larger than the planetary signal at the output 
of the interferometer. Fortunately, the high angular resolution 
provided by the long imaging baseline of Darwin/TPF in the 
X-array configuration is sufficient to spatially distinguish most 
of the extended exozodi emission from the planetary signal and 
only the hole in the dust distribution near the planet significantly 
contributes to the noise level. Considering the full wavelength 
range of Darwin/TPF, we show that the tolerable dust density is 
about 15 times the solar zodiacal density for face-on systems and 
decreases with the disc inclination. In practice, this constraint 
might be relaxed since we examined a resonant ring model that 
does not include dust from highly eccentric or inclined parent 
bodies, the effects of grain-grain collisions, or perturbations by 
additional planets, all of which can reduce the contrast of the 
resonant ring and improve the tolerance to the exozodiacal dust 
density. 

These results show that asymmetric structures in exozodi- 
acal discs around nearby main sequence stars are one of the 
main noise sources for future exo-Earth characterization mis- 
sions. A first solution to get around this issue is to have a long 
imaging baseline architecture which resolves out the more spa- 
tially extended emission of the exozodiacal cloud from the point- 
like emission of planets. The stretched X-array configuration is 
particularly convenient in that respect. The second solution is 
to observe in advance the nearby main sequence stars and re- 
move from the Darwin/TPF target list those presenting a too 
high dust density or disc inclination. The FKSI nulling inter- 
ferometer would be ideal in that respect with the possibility to 
detect exozodiacal discs down to the density of the solar zo- 
diacal cloud. Ground-based nulling instruments like LBTI and 
ALADDIN would also be particularly valuable. 

Appendix A: Deriving instability noise constraints 

We address here the instability noise and derive the constraints 
on the instrument stability required to detect an Earth-like planet 
orbiting at 1 AU around a G2V star loca t ed at 15 pc. The analy- 
sis follows the analytical method of lLavl d2004l) which was orig- 
inally applied to the DCB configuration at 10/i/m. The goal is to 
extend the study to the X-array architecture and to short wave- 
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Fig. A.l. Instability noise with respect to the rms OPD errors for 
different array architectures assuming 1/f-type PSDs and a rms 
amplitude mismatch of 0. 1 % (defined on the 0- 1 4 Hz frequency 
range). The level of shot noise is represented by dotted curves for 
each configuration. The figure has been plotted for 4-m aperture 
telescopes operating at lOyum and a G2V star located at 15 pc 
(surrounded by an exozodiacal cloud of 1 zodi). 



lengths where instability noise is the most dominant. We define 
the limiting rms OPD and amplitude errors such that instability 
noise is dominated by a factor 5 by shot noise over a single rota- 
tion of 50000 s. Assuming that instability noise is totally uncor- 
rected with the rotation angle, this factor stays unchanged over 
multiple rotations. In practice, instability noise can be correlated 
with the rotation angle due to perturbations such as solar heating 
effects but it should be possible to remove them by measuring 
and correcting the amplitudes and phases at intervals during the 
rotation. 

For sake of comparison with 03(2004), we consider first 
4-m aperture telescopes operating at lOyt/m and extend the re- 
sults to 2-m aperture telescopes operating at 7 pm in Table I A. fl 
Assuming a spectral resolution of 20, Fig. I A. 1 1 shows instabil- 
ity noise with respect to the rms OPD error for three different 
configurations: the DCB as defined in lLavl (12004 . the X- array 
with a 1:2 aspect ratio (X-ar2) and the X-array with a 1:6 as- 
pect ratio (X-ar6). These curves have been computed assuming 
1/f-type PSDs for amplitude and OPD errors with rms values 
defined on a frequency range from l/t lot to 10 4 Hz, where t rot is 
the rotation period, and a G2V star located at 15 pc. Shot noise 
is represented by dotted curves. It is higher for the DCB due to 
geometrical leakage (the nulling baseline is twice larger than for 
the X-array) and presents a slight increase for large rms OPD er- 
rors due to instrumental leakage. Instability noise is also higher 
for the DCB than for the X-array. This is because the planetary 
signal is mostly modulated at lower frequencies where the in- 
stability noise is higher for 1/f-type noises. This is illustrated in 
Fig. |A. 21 showing the chopped planet detection rate with respect 
to the rotation angle of the array (upper figure) and the corre- 
sponding Fourier amplitudes (lower figure). 

Fig. I A. 1 1 shows that shot noise dominates instability noise 
by a factor 5 for rms OPD errors of about 0.5 nm, 0.5 nm and 
0.6 nm respectively for the DCB, X-a r2 and X-ar 6. The slight 
discrepancy that can be mentioned with Lav (2004) is due to two 



; b. 






■ : X-ar2 ; 

■ : X-or6 - 

■ : DCB - 


Li 


I 




Ji.. 


L i 


1 , 



5 10 15 20 

Harmonic of rotation frequency 



Fig. A.2. Upper: Chopped planet detected photon rate as a func- 
tion of array rotation angle for the different architectures. The 
planet is assumed to be located at 47 mas from a G2V star lo- 
cated at 15 pc. Lower: Corresponding Fourier amplitudes. Only 
odd harmonics are present because of phase chopping. 



factors. In addition to the instrument throughput of 10% used in 
EiyJ([2004), our study accounts for the coupling efficiency (about 
72% for the on-axis light) and for the quantum efficiency of the 
detectors (70 %). We also combine the two chop states whereas 
the results of lLavl d2004l) are given for only one. Although in- 
stability noise is higher for the DCB configuration than for the 
X-array, the limiting rms opd error is of the same order due to 
the higher shot noise. Because instability noise is directly pro- 
portional to the stellar flux, the constraints are even more strin- 
gent at 7-pm where the star is brighter and the level of shot noise 
basically the same than at lOyum (see Table lA.JV Considering an 
rms amplitude mismatch of 0.1%, the OPD has to be controlled 
to a level of 0.3 nm rms for the three considered configurations. 
These constraints are slightly relaxed for 2-m aperture telescopes 
because shot noise is relatively more dominant (shot noise is pro- 
portional to the squared root of the stellar flux while instability is 
directly proportional to the stellar flux). For the results at 7 pm, 
we use the same baseline length of 20 m with the same rotation 
period of 50000 s. These requirements are very stringent and are 
only marginally compliant with state-of-the-art active control, so 
that potential ways to mitigate the harmful effect of instability 
noise have been investigated. 

- A first solution, proposed by lLavl d2006f) . consists in stretch- 
ing the array and applying a low-order polynomial fit to the 
instability noise (as a function of wavelength). By stretch- 
ing the array, i.e. increasing the imaging baseline of the 
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Table A.l. Limiting rms OPD values computed for a G2V star located at 15 pc such that shot noise dominates instability noise by a 
factor 5. The results are given at 7 fim and 10/im, assuming a spectral resolution of 20, for different rms amplitude mismatches and 
shapes of the instability noise power spectrum. Three configurations are shown, the X-array with an a spect ratio of 2 (X-ar2), the 
X-array with an aspect ratio of 6 (X-ar6) and the Dual Chopped Bracewell (DCB, configuration used in Lav 2004b. 
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X-array, the interference pattern orthogonal to the nulling 
pattern shrinks. As the modulation map scales with wave- 
length, the planetary signal transmitted by the interferome- 
ter will then be a rapidly oscillation function of wavelength. 
On the other hand, instability noise is shown to be a low- 
order polynomial of the optical frequency (1/A). Therefore, 
by removing a low-order polynomial fit to the detected sig- 
nal as a function of wavelength, the instability noise contri- 
bution is efficiently subtracted while preserving most of the 
planetary signal. This operation can actually be performed 
directly in the cross-correlation, by using a modified planet 
template where the polynomial components have been re- 
moved. Because this method strongly relies on the separa- 
tion of the nulling and imaging baselines, it can only be effi- 
ciently applied with the X-array architecture. 

- Another solution, based on th e coherence prop erties of 
starlight, has been proposed by lLane et ail d2006l) to sepa- 
rate the contributions from the planet and the instrumental 
leakage. The idea is to mix the electric fields of the leakage 
with that of a separate reference beam, also from the star, in 
order to form fringes (as long as the relative path delays are 
maintained within the coherence length). The light from the 
planet, being not coherent with the starlight, will not form 
fringes. Using as reference beam the bright output of a pair- 
wise nulling beam combiner, the amplitude and phase mis- 
matches in the input beams can be extracted from the fringe 
pattern, allowing the reconstruction of the the stellar leakage. 

- A third solution has recently been proposed and tested at IAS 
(Gab or et alj |2008). The principle is to successively apply 
two (< A/2) opposite OPD offsets to one of the beams in or- 
der to derive the actual position of the minimum in the trans- 
mission map. In this way, one feeds back the actual OPD 
errors to the delay line and prevents drifts from appearing. 
The same principle can be used to avoid amplitude drifts, by 
either blocking all but one beam to measure its actual am- 
plitude or by modulating its amplitude as in the case of the 
OPD. The frequency at which this process is carried out de- 
pends on the input power spectra of the noise sources. It has 
been demonstrated experimentally in the lab that this process 
efficiently suppresses the 1/f-type noise generally present in 
the stellar contribution at the output of a Bracewell inter- 
ferometer. This technique can theoretically be applied to any 
nulling configuration, but its efficiency decreases as the num- 
ber of beams increases. 



Others possible approaches (e.g. application of closure phase 
techniques, new chopping processes, exploiting correlations be- 
tween measurements taken at different wavelength bins) have 
been suggested but still need to be investigated. However, these 
very stringent requirements have been derived with the assump- 
tion that the instability noise present an 1/f-type power spectrum. 
In fact, the shape of the instability noise spectrum strongly de- 
pends on the input spectrum of noise fluctuations. Pure 1/f noise 
might in fact turn out to be a very pessimistic assumption. Given 
that the Darwin/TPF system is designed to have three levels of 
control loops that manage the OPD, and two control loop lev- 
els for tip/tilt, the resulting residual phase and intensity fluc- 
tuations that are responsible for the instability noise are most 
likely to have a flat input spectrum. This has a strong influence 
on the magnitude of the phase-amplitude cross terms, which are 
the main contributor to instability noise. For instance, the tol- 
erable rms OPD errors for a white spectrum defined on the [0- 
10 Hz] frequency range are significantly relaxed as indicated in 
Table I A. ll In particular at low frequencies, where the planetary 
signal is modulated and the instability noise arises, the 1/f spec- 
trum diverges, while the flat spectrum increases linearly. In prac- 
tice, the situation might even be better, as predictive or Kalman 
filtering may be used to remove low-frequency power. A compli- 
cation, in particular for the micro Hz domain, is that zero point 
drifts of the control loops, e.g. related to electronic drifts in the 
sensors, remain undetected, and thus uncorrected. They will re- 
sult in a low-frequency 1/f component. A simple solution is to 
switch the incoming beams with respect to the control loops. 
The reason why this could work is that the extremely high pre- 
cision required for Darwin/TPF is only on relative quantities, 
namely, the phase and amplitude differences between any two 
beams. Absolute offsets of phase and amplitude, that apply to 
all three beams simultaneously, do not contribute to the instabil- 
ity noise, since they do not affect the null. The beams must be 
switched with respect to the control inputs with a time constant 
over which the drift can be considered constant. In summary, 
even though at present a post-processing method has not been 
identified, it is highly likely that clever and careful engineering 
may already reduce instability noise to harmless levels. 

Acknowledgements. The authors are grateful to Lisa Kaltenegger (Harvard- 
Smithsonian Center for Astrophysics) for providing the updated Darwin 
catalogue, Jean Surdej (IAGL), Arnaud Magette (IAGL), Dimitri Mawet 
(NASA/JPL), Peter Lawson (NASA/JPL), Oliver Lay (NASA/JPL), Pierre Riaud 
(IAGL) and Virginie Chantry (IAGL). This research was supported by the 
International Space Science Institute (ISSI) in Bern, Switzerland ("Exozodiacal 
Dust discs and Darwin" working group, http://www.issibern.ch/teams/exodiist/). 



D. Defrere et al.: Impact of exozodiacal clouds on the performance of future life-finding space missions 



17 



DD and CH acknowledge the support of the Belgian National Science 
Foundation ("FRIA"). OA acknowledges the support from a F.R.S.-FNRS 
Postdoctoral Fellowship. DD and OA acknowledge support from the 
Communaute francaise de Belgique - Actions de recherche concertees - 
Academie universitaire Wallonie-Europe. 



References 

Absil, O. 2001, Master's thesis, Liege University, Liege, Belgium, 

http://www. aeos. ulg. ac.be/ 
Absil, O. 2006, PhD thesis, Liege University, Liege, Belgium 
Absil, O., Defrere, D., Coude du Foresto, V., et al. 2008a, in Proc. SPIE, Vol. 

7013 

Absil, O., di Folco, E., Merand, A., et al. 2006, A&A, 452, 237 
Absil, O., di Folco, E., Merand, A., et al. 2008b, A&A, 487, 1041 
Absil, O., Mennesson, B., Le Bouquin, J.-B., et al. 2009, ArXiv e-prints 
Akeson, R. L., Ciardi, D. R., Millan-Gabet, R., et al. 2009, ApJ, 691, 1896 
Alonso, R., Auvergne, M., Baglin, A., et al. 2008, A&A, 482 
Angel, J. R. P. & Woolf, N. J. 1997, ApJ, 475, 373 
Barge, P., Baglin, A., Auvergne, M., et al. 2008, A&A, 482 
Beichman, C. A., Bryden, G., Stapelfeldt, K. R., et al. 2006, ApJ, 652, 1674 
Borucki, W. J., Koch, D. G., Lissauer, J., et al. 2007, in Astronomical Society 
of the Pacific Conference Series, Vol. 366, Transiting Extrapolar Planets 
Workshop, ed. C. Afonso, D. Weldrake, & T. Henning, 309-+ 
Bracewell, R. N. 1978, Nature, 274, 780 

Carle, E. 2005, EMMA configuration: evaluation of optical performances, Tech. 

Rep. Issue 1, ESA (SCI- A/279) 
Chazelas, B., Brachet, F, Borde, P., et al. 2006, Appl. Opt., 45, 984 
d'Arcio, L. 2005 

Defrere, D., Absil, O., Coude Du Foresto, V., Danchi, W. C, & den Hartog, R. 

2008a, A&A, 490, 435 
Defrere, D., Lay, O., den Hartog, R., & Absil, O. 2008b, in Proc. SPIE, Vol. 7013 
den Hartog, R. 2005a, DARWIN science performance prediction, Tech. Rep. 

Issue 1, ESA (SCI- A/300) 
den Hartog, R. 2005b, The DARWINsim science simulator, Tech. Rep. Issue 1 , 

ESA (SCI-A/297) 

Dermott, S. E, Jayaraman, S., Xu, Y. L., Gustafson, B. A. S., & Liou, J. C. 1994, 
in, 7 19-723 

Dermott, S. F, Nicholson, P. D., Burns, J. A., & Houck, J. R. 1985, in 
Astrophysics and Space Science Library, Vol. 119, IAU Colloq. 85: Properties 
and Interactions of Interplanetary Dust, ed. R. H. Giese & P. Lamy, 395—409 

Di Folco, E., Absil, O., Augereau, J.-C, et al. 2007, A&A, 475, 243 

Dohnanyi, J. S. 1969, J. Geophys. Res., 74, 2531 

Fridlund, C. V. M., d'Arcio, L., den Hartog, R., & Karlsson, A. 2006, in Proc. 
SPIE, Vol. 6268 

Fridlund, M. 2005, Darwin Science Requirements Document, Tech. Rep. Issue 

5, ESA (SCI-A) 
Gabor, P., Chazelas, B., Brachet, F, et al. 2008, A&A, 483 
Gillon, M., Pont, F, Demory, B.-O., et al. 2007, A&A, 472, L13 
Greaves, J. S., Holland, W. S., Wyatt, M. C, et al. 2005, ApJ, 619, L187 
Hinkley, S., Oppenheimer, B. R., Soummer, R., et al. 2007, ApJ, 654, 633 
Hinz, P. M., Bippert-Plymate, T., Breuninger, A., et al. 2008, in Proc. SPIE, Vol. 

7013 

Kalas, P., Graham, J. R., & Clampin, M. 2005, in , 1067-1070 
Kaltenegger, L., Eiroa, C, & Fridlund, C. V. M. 2008, ArXiv e-prints 
Karlsson, A. L., Wallner, O., Perdigues Armengol, J. M., & Absil, O. 2004, in 

Proc. SPIE, ed. W. A. Traub, Vol. 5491, 831-+ 
Kelsall, T, Weiland, J. L., Franz, B. A., et al. 1998, ApJ, 508, 44 
Ksendzov, A., Lay, O., Martin, S., et al. 2007, Appl. Opt., 46, 7957 
Ksendzov, A., Lay, O., Martin, S., et al. 2008, Appl. Opt., 46, 7957 
Kuchner, M. J. & Holman, M. J. 2003, ApJ, 588, 1110 
Landgraf, M. & Jehn, R. 2001, Ap&SS, 278, 357 
Lane, B. F, Muterspaugh, M. W., & Shao, M. 2006, ApJ, 648, 1276 
Lawson, P. R., Lay, O. P., Martin, S. R., et al. 2008, in Proc. SPIE, Vol. 7013 
Lay, O. P. 2004, Appl. Opt., 43, 6100 
Lay, O. P. 2005, Appl. Opt., 44, 5859 
Lay, O. P. 2006, in Proc. SPIE, Vol. 6268 

Lay, O. P., Martin, S. R., & L., H. S. 2007, in Proc. SPIE, Vol. 6693 
Leger, A. & Herbst, T. 2007, ArXiv e-prints, 707 

Leger, A., Rouan, D., Schneider, J., Alonso, B., & Samuel, E. subm., A&A 
Liou, J.-C. & Zook, H. A. 1999, AJ, 1 18, 580 

Macintosh, B. A., Becklin, E. E., Kaisler, D., Konopacky, Q., & Zuckerman, B. 

2003, ApJ, 594, 538 
Mawet, D., Hanot, C, Lenaers, C, et al. 2007, Optics Express, 15, 12850 
Mayor, M., Bonfils, X., Forveille, T., et al. 2009, A&A, in press 
Mayor, M. & Queloz, D. 1995, Nature, 378, 355 
Mennesson, B. & Mariotti, J. M. 1997, Icarus, 128, 202 



Moro-Martin, A. & Malhotra, R. 2005, ApJ, 633, 1150 
Moutou, C, Pont, F, Barge, P., et al. 2005, A&A, 437, 355 
Ozernoy, L. M., Gorkavyi, N. N, Mather, J. C, & Taidakova, T. A. 2000, ApJ, 
537, L147 

Peters., R., Gappinger, R., Lawson, P., & Lay, O. 2009, JPL, Document D-60326 
Reach, W. T, Franz, B. A., Weiland, J. L., et al. 1995, Nature, 374, 521 
Reche, R., Beust, H., Augereau, J.-C, & Absil, O. 2008, A&A, 480, 551 
Richardson, L. J., Deming, D., Horning, K., Seager, S., & Harrington, J. 2007, 

Nature, 445, 892 
Robertson, H. P. 1937, MNRAS, 97, 423 

Roques, F, Scholl, H., Sicardy, B., & Smith, B. A. 1994, Icarus, 108, 37 
Ruilier, C. & Cassaing, F. 2001, Journal of the Optical Society of America A, 
18, 143 

Schneider, G., Weinberger, A. J., Becklin, E. E., Debes, J. H., & Smith, B. A. 

2009, AJ, 137, 53 
Serabyn, E. 2009, ApJ, 697, 1334 

Serabyn, E. & Colavita, M. M. 2001, Appl. Opt., 40, 1668 

Stark, C. C. & Kuchner, M. J. 2008, ApJ, 686, 637 

Stark, C. C. & Kuchner, M. J. 2009, ArXiv e-prints 

Stark, C. C, Kuchner, M. J., Traub, W. A., et al. 2009, 41, 501 

Thiebaut, E. & Mugnier, L. 2006, in IAU Colloq. 200: Direct Imaging of 

Exoplanets: Science and Techniques, 547-552 
Wilner, D. J., Holman, M. J., Kuchner, M. J., & Ho, P. T. P. 2002, ApJ, 569, 

L115 



